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EXTERNAL  NONLINEAR  FOCUSING  AND  ITS 
APPLICATION  TO  ULTRAFAST  OPTICAL  POWER  LIMITING 


1 .  INTRODUCTION 


If  a  laser  beam  of  sufficient  power  passes  through  a  thin 
nonlinear  medium  it  will  suffer  an  optical  phase  shift.  In  particular,  when 
a  laser  beam  of  gausslan  cross-section  interacts  with  a  refractive 
nonlinearity,  the  resultant  phase  shift  is  found  fl-3)  to  generate  strong 
interference  effects  in  the  external  far  field.  The  transverse  profile  in 
the  far-field  then  consists  of  a  set  of  concentric  rings  whose  number 
increases  with  beam  power  (Miller  et  al  ill,  Weaire  et  al  12],  Miller  et  al 
(31).  The  present  report  studies  the  power-limiting  capability  of  a  thin 
nonlinear  medium  in  conjunction  with  an  'on-axi3'  far-field  aperture.  The 
same  solutions  apply  if  a  lens  placed  after  the  nonlinear  medium  focuses  the 
phase-shifted  beam  onto  an  aperture  near  its  focal  plane  (see  figure  la).  At 
low  beam  powers  virtually  all  of  the  energy  is  transmitted  through  the 
aperture,  while  at  higher  powers  a  large  part  of  the  energy  is  absorbed  by  the 
intercepting  screen.  If  the  nonlinear  refractive  coefficient,  n2,  is 
positive,  external  self-focusing  (Kaplan  (41)  occurs  and  the  principal  focal 
spot  is  shifted  towards  the  Incoming  beam,  self-def ocusing  takes  place  if  n2 
is  negative,  and  the  diffraction  focus  is  shifted  to  the  far  side  of  the 
linear  focal  plane. 


A  modification  of  the  above  arrangement  focuses  the  radiation  to  a 
spot  within  the  nonlinear  medium  with  a  second  linear  lens  (depicted  in  figure 
lb),  thereby  providing  maximum  intensity  for  the  generation  of  nonlinear 
processes.  Both  continuous  and  pulsed  radiation  can  be  limited;  in  the 
latter  case  the  external  focal  spot  performs  a  rapid  'whiplash'  motion  towards 
the  nonlinear  medium  as  the  pulse  moves  through  the  medium.  Boggess  et  al 
(5)  have  achieved  power  limiting  of  50  psec  Nd : yag  pulses  at  1.064  pm  with 
3uch  a  system,  using  silicon  as  the  nonlinear  medium.  The  mechanism  in  this 
case  is  lrradiance-dependent  def ocusing.  Soileau  et.  al  16-91  employed  a 

similar  experimental  arrangement  to  achieve  power  limiting  of  1.064  and 
0.532  pm  ultrashort  pulses  using  liquid  CS2  as  the  medium.  In  these 
experiments,  however,  the  nonlinear  processes  were  oily  generated  when  the 
laser  power  P^  exceeded  the  3elf-trapping  threshold  power  for  the  liquid. 


1 


In  another  modification,  which  is  yet  to  be  explored  and  tested 
experimentally,  a  second  aperture  is  placed  at  the  focal  position  in  the 
nonlinear  medium.  This  can  perform  several  functions,  for  example  it  can  act 
as  an  efficient  spatial  filter,  or  can  be  of  use  where  an  input  beam  is  off- 
axis  or  for  an  input  profile  consisting  of  a  large  number  of  small  spots. 


2.  CALCULATIONS 


2.1  General 


The  Important  features  in  the  diffraction  pattern  arising  from  the 
useful  arrangement  in  fig.  lb  are  also  produced  by  the  simpler  system  shown  in 
Fig.  la.  Hence  it  is  instructive  to  study  the  model  of  external  self-focusing 
by  the  simpler  optical  arrangement  depicted  in  fig.  la  (see  references  to  the 
theory  do,  ill).  The  method  of  calculation  is  similar  to  that  employed  for 
aberration  calculations  near  the  focal  plane  of  a  linear  lens  1121.  The 
normalised  intensity  distribution  of  a  beam  of  initial  width  wQ  aberrated  by  a 
thin  nonlinear  layer  (thickness  l)  has  been  shown  (Hermann,  (lol)  to  be 
proportional  to  z-2  in  the  vicinity  of  the  low-power  focus,  with 

I  *°  I 

I  (p;<|i,$)  -  ]  J"qf(x) jq(4ix)  exp  jltpF(x)2-$x2 1^  .x  dxi2  (l) 


where  the  following  definitions  apply: 

(a)  The  focal  length  z  has  been  scaled  according  to  ij-  z/R,,  where 

Rd  ”^2  ko  w<>  13  the  diffraction  length  and  kQ  is  the  wavenumber  of 

the  incident  radiation. 

(b)  The  radial  coordinate  in  the  far  field,  r,  has  been  conveniently 

scaled  as  +-  2r/(w  £). 

o’ 

(c)  The  radial  coordinate  at  the  medium,  r',  has  been  scaled  as 

x  -  r'/Wg,  and  F(x)  is  the  shape  of  the  transverse  profile  of  the 
incident  field  amplitude. 

(d)  The  spatial  coordinate  representing  displacement  from  the  linear 

(low-power)  f^cal  plane  in  the  axial  direction,  y,  has  been  scaled 

as  +  -  y  R_,/z^ . 

<1 

(e)  The  nonlinear  coefficient  ,  which  has  been  defined  in  reference  10, 

can  be  regarded  as  a  parameter  in  the  present  model .  It  can  be  written 
as  p  -  v  where  v  is  a  beam-shape  parameter  which  measures  the 

transverse  spread  of  the  incident  beam  and  Px  is  a  scaling  power  given 
by  <3ee  ref.  10)  P.^  -  q  ji  -  exp( -c^l) kQn2 f ~  . 

Here,  is  the  linear  absorption  coefficient  and  n0,n2  are  the 
linear  and  nonlinear  refractive  indices  respectively.  The 
coefficient  p  is  an  intensity-dependent  quantity. 

(g)  Jgl+x)  is  the  zero-order  Bessel  function. 
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Note  that  the  intensity  of  the  diffracted  beam  has  been  scaled  to  the  maximum 
input  intensity,  so  that  I  remains  finite  even  when  p-  o.  The  displacement 
coordinate  $  has  been  formed  with  due  allowance  (to  lowest  order)  for  the 
effect  of  a  linear  lens  upon  the  phase  of  the  optical  field.  From  Fourier 
optics  it  is  known  that  the  field  is  altered  by  a  phase  factor  exp(ix2/  ). 


The  amplitude  function  F(x)  can  be  expanded  for  a  single  radlaily- 
symmetric  spot  as  a  power  series  of  the  form 

f(x)  -  1  -  x2  +  hx4  +  . . .  (2) 

Substituting  equation  (2)  into  equation  (1)  gives 


Kp;*,4>> 


f(x) 


J0<(|ix)  exp.{l(p(l 


+  2b)x4  -  (2p  +  $)x2)}.x  dx 


(3) 


which  displays  a  power-dependent  displacement  of  the  various  profile  features 
towards  the  focusing  system.  This  type  of  behaviour  represents  an  application 
of  the  'displacement  theorem'  (reference  1 1 2 1 ,  p  462).  Furthermore,  there  is 
a  power-dependent  modification  to  the  phase,  which  is  at  least  of  the  fourth 
degree  in  x,  producing  in  turn  a  power-dependent  aberration  of  the  various 
profile  features. 


The  sign  of  determines  whether  the  system  is  self-focusing 
(positive)  or  self-defocusing  (negative).  Equation  (l)  reveals  that  changing 
the  sign  of  p  is  equivalent  to  reflecting  the  intensity  distribution  about  the 
linear  focal  plane,  le., 


l(-p;<|i,<fr)  -  l( p;i|i, -$)  (4) 

The  defocusing  calculation  is  therefore  equivalent  to  the  self-focusing 
calculation,  if  the  spatial  domain  is  restricted  to  the  region  of  the  linear 
focal  plane. 


2 . 2  Single-mode  input 


We  have  Integrated  equation  (l)  numerically  for  different  laser 
input  profiles  and  have  plotted  the  normalised  intensity  distribution  in 
the  region  of  the  linear  focal  plane  as  a  function  of  The  integration 

wa3  performed  with  the  U3e  of  adaptive  routines  which  used  both  Gauss  10-point 
and  Konrod  21-polnt  rules  (131.  The  numerical  accuracy  of  the  integration 
technique  was  predetermined  to  be  10-5,  a  level  considered  to  be  satisfactory, 
and  as  an  additional  check  on  the  overall  accuracy,  the  transmitted  laser 
power  for  a  nonabsorbing  system  was  calculated  in  the  linear  focal  plane  and 
was  found  to  be  conserved  to  better  than  one  percent.  The  results  are 
depicted  in  figures  2  and  3  for  single-spot  profiles,  with  gaussian  and 
parabolic  shapes  respectively,  and  with  increasing  positive  values  of  p.  in 
the  gaussian  case,  the  shape  of  the  focal  spot  in  the  limit  p-»-o  is  easily 
calculated  as 

I(0;ij>,4»  -  |(1  +  4,2)-1  exp  {-  J  i|p2/(i  +  (5) 
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which  displays  a  gaussian  shape  transversely  and  a  lorentzlan  shape  axially. 
For  increasing  p  the  focal  region  stretches  out  along  the  optic  axis,  and  its 
centre  moves  away  from  the  linear  focal  plane.  Furthermore,  side  peak3 
develop  and  grow  as  p  increases,  and  eventually  the  focal  ridge  itself 
develops  into  a  sequence  of  peaks,  or  bright  spots.  It  is  possible  to 
evaluate  the  diffraction  integral  on-axis  in  terms  of  a  confluent 
hypergeometric  function  (which  also  can  be  represented  in  terms  of  generalised 
Fresnel  integrals )  ( i 4 ) : 


Kp;o,+)  -  jd  +  -fr2)-1 


jFjd,  |  (14>+  3);  -  ip) 


(6) 


Note  that  the  property  given  by  eqn  U)  is  also  implicit  in  eqn  (6).  This  on- 
axis  formula  possesses  the  property  that  the  lorentzian  shape  becomes  multiply 
peaked  as  p  increases  from  zero.  Similar  structures  associated  with  len3 
aberrations  (15)  have  been  held  responsible  for  effects  produced  during  laser- 
induced  breakdown  in  gases  (IS),  such  as  separated  regions  of  Intense 
ionization  along  the  optic  axis  near  the  lens  focus. 


It  is  also  apparent  from  figures  2a-c  that  the  intensity  along  the 
optic  axis  declines  rapidly  beyond  the  linear  focal  plane,  suggesting  that  for 
the  purpose  of  maximising  the  power-limiting  capability  of  a  real  self- 
focusing  device  it  is  desirable  to  position  an  aperture  slightly  beyond  the 
linear  focal  plane  (or  before  that  plane  for  a  defocusing  device).  In  the 
experiments  reported  by  reference  3  it  was  found  that  the  largest  intensity 
peaks,  in  the  far-fleld  profile  generated  by  a  gaussian-like  input  with  a 
defocusing  medium,  can  reside  in  the  wings  of  the  profile  when  p  is  large 
enough.  This  behaviour  has  been  observed  also  in  the  present  calculations  if 
the  plane  of  observation  13  positioned  before  the  linear  focal  plane.  At  the 
focal  plane,  the  insertion  of  a  precisely  gaussian  input  profile  into  the 
diffraction  integral  leads  to  diffraction  patterns  in  which  the  largest  peak 
is  always  on-axis  (101.  This  suggests  that  the  input  profiles  claimed  in 
reference  3  might  not  have  been  truly  gaussian,  but  truncated  gausslans,  for 
which  the  higher-order  terms  are  different.  Noteworthy  in  thi3  connexion  is 
the  observation  that  incident  laser  beams  possessing  field  amplitudes  of 
parabolic  cross-section  behave  somewhat  differently  in  the  present  model  to 
those  with  precisely  gaussian  profiles;  in  the  parabolic  case  the  features 
along  the  optic  axis  all  translate  together  with  respect  to  increasing  p. 
Hence,  for  large  enough  p,  any  plane  cutting  the  optic  axis  perpendicularly  in 
the  neighbourhood  of  the  linear  focal  plane  intercepts  an  Intensity  profile 
with  peak  maxima  in  the  wings,  rather  than  on-axis.  Further  experimental 
work  is  evidently  desirable  so  that  these  two  explanations  for  the  shift  of 
power  away  from  the  axis  to  the  wings  can  be  satisfactorily  disentangled,  and 
also  to  establish  the  relevance  of  higher-order  effects  (represented  by  the 
wings  of  the  beam)  to  the  qualitative  features  of  the  diffraction  pattern. 


2.3  Higher-mode  input 


More  compiicated  behaviour  has  been  found  <ith  models  in  which  the 
laser  input  profile  takes  the  form  of  a  symmetric  spot  surrounded  by  one  or 
more  concentric  rings,  or  of  concentric  rings  without  a  central  spot.  The 
annular  input  profile  given  by  F(x)  -  4x2(i-x2)  with  o  4x  4  l  has  been 
explored,  and  the  focal  distributions  are  shown  in  figures  4a,b,c.  in  this 
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case  there  is  a  focal  distribution  on  both  sides  of  the  linear  focal  plane; 
as  p  increases  the  focal  region  becomes  elongated  and  attenuated,  and  internal 
peaks  form.  The  number  of  peaks  in  the  focal  region  increases,  and  the  peaks 
move  apart,  as  p  increases.  This  combination  of  both  self-focusing  and  self- 
defocusing  occurs  because  the  velocity  of  the  electromagnetic  wave  within  the 
nonlinear  medium  is  diminished  in  regions  of  high  intensity,  so  that  a  convex 
intensity  profile  renders  the  beam  convergent  while  a  concave  profile  produces 
divergent  behaviour.  Thus,  an  annular  input  beam  traversing  a  nonlinear 
medium  generates  a  wavefront  with  both  converging  and  diverging 
characteristics,  giving  large  peaks  on  each  side  of  a  linear  focal  plane. 

Note  that  amplitude  profiles  of  the  form  F(x)  -  ax2  -  bx4  +...  possess  a 
symmetric  distribution  along  the  optic  axis  with  respect  to  the  origin  at  the 
focal  plane.  The  reason  is  evident  from  an  analysis  analogous  to  that  used  to 
obtain  equation  (3).  Asymmetries  appearing  off-axis  are  attributable  to  terms 
of  order  x4  and  higher.  Annular  amplitudes  which  are  odd  in  x  will  generally 
suffer  displacement  of  the  focal  region  (in  proportion  to  p)  as  well  as 
asymmetries  along  the  optic  axis. 


2.4  Inclusion  of  nonlinear  absorption 

The  extension  of  the  above  calculations  to  incorporate  nonlinear 
absorption  (below  saturation)  has  been  carried  out,  following  the  theoretical 
analysis  of  reference  10.  Instead  of  equation  (l),  we  now  evaluate 


I(a, 


F(x)J  (  x)  . 

- 9— — !/  exp.  ji(ain(l+pF(xr>  -<t>x2]} 

{l+  F<x)2}  /2 


x  dx 


(7) 


where  a  is  essentially  the  ratio  of  nonlinear  refraction  to  nonlinear 
absorption  coefficients,  and  p  is  the  power  parameter  appropriately  scaled. 
Note  that  p  is  also  proportional  to  the  nonlinear  absorption  coefficient. 
Changing  the  sign  of  n2  (and  therefore  a)  is  again  equivalent  to  reflecting 
the  intensity  distribution  given  by  equation  (7)  about  the  linear  focal 
plane. 


we  have  calculated  the  focal  intensity  distributions  for  a  parabolic 
input  (F(x)  -  1  -  x2)  with  a  fixed  refraction  parameter  and  several  nonlinear 
absorption  coefficients.  The  main  feature  to  emerge  from  these  calculations 
is  the  opposing  behaviour  of  the  absorptive  and  refractive  mechanisms  towards 
focal  displacement.  Where  nonlinear  absorption  dominates,  the  diffraction 
focus  remains  at  the  linear  focal  plane,  suitably  attenuated  in  response  to 
the  magnitude  of  the  nonlinear  absorption  coefficient.  Where  nonlinear 
refraction  dominates,  the  diffraction  focus  is  displaced  from  this  plane  to 
the  maximum  extent.  This  interplay  between  the  two  mechanisms  explains  many 
of  the  power-limiting  features  reported  in  reference  2  and  the  experimental 
results  for  the  silicon  switch  reported  in  references  (5). 
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3. 


DISCUSSION 


The  present  model  embraces  both  nonlinear  refraction  and  absorption 
within  the  active  medium.  It  adopts,  for  an  effective  dielectric 
susceptibility  X  of  the  medium,  the  expression 


X(|e|2)  ocXg  +  iQj  +  (kQn2/nQ  *  ia2*  lE|2  (8) 

where  E  Is  the  envelope  of  the  input  field  amplitude,  and  a is  the  nonlinear 
absorption  coefficient.  Higher-order  terms,  and  hence  saturation  effects, 
have  been  neglected.  The  model  is  generally  Inadequate  for  describing  the 
interaction  of  pulsed  radiation  with  materials  whose  refractive  indices  are 
strongly  Influenced  by  photogenerated  free  carriers  (assuming,  of  course,  that 
the  pulse  width  does  not  exceed  carrier  lifetimes).  Dynamic  features  have  to 
be  allowed  for  in  these  cases,  and  time-dependence  must  be  admitted  to  the 
theoretical  description  both  within  the  wave  equation  and  through  the  time 
evolution  of  the  free-carrier  populations.  Despite  the  absence  of  these 
features  from  the  model,  there  is  agreement  for  appropriate  values  of  $  and  p 
with  the  features  reported  in  previous  experimental  work  for  both  continuous 
radiation  (middle  IR)  and  pulsed  radiation  (near  IR) .  Theoretical  work 
giving  more  detailed  quantitative  predictions  will  await  the  development  of 
satisfactory  dynamic  models  for  absorption  and  refraction  (a)  by  laser¬ 
generated  plasmas  within  liquids,  and  (b)  for  multiphoton  and  free-carrier 
effects  within  solid  media. 


Saturation  (or  partial  saturation)  of  the  nonlinear  absorption  is 
likely  at  higher  laser  powers.  Recent  research  by  Hill  et  al  (16)  with  a 
gaussian  beam  at  wavelength  10.6  pm  traversing  a  crystal  of  cadmium  mercury 
telluride  at  ip0Kj^as  shown  a  nonlinear  refractive  index  change 

The 


-yj  x  10  "  I  where  I oc  |  E|‘  is  the  incident  intensity  in  w/cm  . 

dependence  reflects  the  important  role  of  Auger  recombination  (see 


references  in  (161)  in  limiting  the  saturation  of  absorption  at  the  band  edge, 
and  similar  behaviour  can  be  expected  with  other  materials  such  as  indium 
antimonide.  The  optical  model  that  has  been  explored  by  these  authors 
differs  from  ours  in  that  they  constrain  the  nonlinear  medium  with  a  circular 
aperture,  and  also  by  their  imposition  of  a  quadratic  approximation  in  the 
description  of  the  phase  profile.  The  latter  simplif ication  has  been 
justified  by  them  on  the  basis  that  saturation  effects  broaden  the  phase 
profile  relative  to  the  intensity  profile.  It  is  of  interest  to  extend  eqn. 
(6)  for  the  gau3sian-input  profile  and  an  infinite  aperture,  by  suitably 
incorporating  a  n^w  parameter  p,  to  account  for  saturation  effects. 

Writing  eyp(-2x  /p  )  for  the  phase  profile,  we  find  instead  of  eqn.  (6): 


I(p;0,$)  -  i  (1  +  $2)  1  j  lFl(l,  1  +  |  p2(l  +  i$);  -ip) 


(9) 


The  corresponding  result  restricted  by  the  phase-profile  approximation 
p(i  -  2x  /p  )  and  by  a  finite  circular  aperture  of  radius  R(<p)  is 
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i(p;o,<fr) 


J*R  exp.  |-x2(l  +  i($  +  2p/H2) )  (  .  x  dx 


11  +  ($  +  2p/p2)l  1  .  e  R  .  islnh2(-^R2)  +  sin2  -jR2($  +  2p/p2)^ 


(10) 


The  displacement  of  the  diffraction  focus  (the  tallest  peak)  is  given  by 


*max  '  -  ' 


a3  anticipated,  and  the  remaining  peaks  within  the  focal  region  are 
symmetrically  disposed.  The  Intensity  distribution,  given  by  equation  (9), 
exhibits  pronounced  asymmetry  along  the  optic  axis  and  is  clearly  shown  in  the 
numerical  results  in  figs.  2a-c.  Furthermore,  eqns.  (9)  and  (10)  differ  in 
the  limit  p— *-0,  except  when  R— *■  .  We  conclude  that  great  care  should  be 
taken  with  the  mathematical  modelling  to  ensure  that  effects  associated  with 
the  tails  of  the  beam  and  phase  profiles  are  fully  accounted  for. 


The  range  of  validity  in  employing  a  quadratic  approximation  for  the 
phase  profile  also  can  be  tested  by  examining  the  case  of  parabolic  input 
amplitude.  Approximating  the  phase  p(l-x2>  ,  where  X-  2/p2,  by  the  first- 
order  expression  p(l-2x2/p2),  we  have 

I  (  p ;  0 , 4> )  =  jj"1  (l  -x2)  exp.  |  -i0x2  }  x  dx|  2 

-  ^  9-4  |  1  -  cos  9-0  sin  0  +  ^  02  ( 11 ) 

2 

where  0  -  $  +  2p/p  .  once  again,  there  is  symmetrical  distribution  about  the 
peak  at  0-0.  The  restriction  to  the  lowest  order*  effect  (displacement  of 
the  focal  spot)  is  clearly  a  more  serious  source  of  error  for  input  profiles 
possessing  substantial  tails. 


In  reference  10  the  power-transfer  curve  (relating  incident  to 
transmitted  power  at  the  steady  state)  was  calculated  for  a  parabolic  input 
profile.  The  latter  shape  was  chosen  by  reason  of  the  relative  computational 
ease,  compared  to  the  gaussian.  There  is  a  large  initial  'hump'  on  the 
parabolic  transfer  curve,  whereas  the  corresponding  gaussian  curve  exhibits  a 
much  less  pronounced  hump.  This  difference  in  behaviour  reflects  the 
different  responses  of  the  respective  focal  peaks  towards  increasing  p. 

In  order  to  Inspect  the  origins  of  the  difference  most  easily,  we  have  drawn 
lsophote3  (lines  of  constant  Intensity)  for  the  two  cases  in  figs.  5-6.  The 
displacement  and  aberration  effects  are  clearly  displayed.  Cross-sections 


The  term  'order'  here  refers  to  powers  of  x2  in  the  series  expansion  for 
the  optical  phase  profile,  and  should  not  be  confused  with  the  more 
general  meaning  of  'order'  in  optical  aberration  theory. 
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perpendicular  to  the  optic  axis  are  shown  for  the  gaussian  case  in  fig.  7,  and 
reveal  a  substantial  relative  redistribution  of  radiant  power  away  from  the 
region  of  the  optic  axis  as  the  intercepting  plane  crosses  the  linear  focal 
plane.  Another  noticeable  difference  in  behaviour  between  the  two  cases  is 
that  the  focal  ridge  for  the  gaussian  input  remains  'anchored',  with  one  end 
always  extending  as  far  as  the  linear  focal  plane,  whereas  for  parabolic  input 
all  of  the  axial  features  translate  together  as  p  increases. 


We  have  been  able  to  effect  a  semi -quantitative  comparison  between 
theory  and  experiment,  despite  the  absence  of  scaling  information.  The 
parameters  relating  to  the  'far-field'  experiment  with  indium  antimonlde 
reported  in  references  1-2  were  inserted  into  our  computer  programme  for 
evaluating  eqn.  (l),  and  intensity  profiles  were  generated  near  the  linear 
focal  plane  for  various  values  of  $.  It  was  found  that  the  best  fit  could  be 
obtained  for  three  different  values  of  p  (specifically  p  =  2.20,  5.32,  8.72) 
if  the  observational  plane  intercepts  the  optic  axis  at  $  =  0.6,  and  a  typical 
result  is  shown  in  fig.  8.  Issues  relating  to  the  shifting  of  radiant  power 
away  from  the  axial  region  should  be  unambiguously  resolved  by  carrying  out 
further  experiments.  For  example,  an  experiment  with  an  Indium  antimonide 
layer  and  material  parameters  corresponding  to  p  -  10  appears  feasible;  the 
extent  to  which  the  present  model  with  gaussian  input  (fig.  2b)  describes  this 
situation  could  be  tested  by  performing  spatial  scans  at  various  intercepting 
planes.  Input  profiles  could  also  be  varied,  using  transversely-graded 
filters  (apodlzatlon) . 

The  various  types  of  behaviour  described  in  section  2  are  highly 
relevant  to  the  performance  and  design  characteristics  of  optical  switches  and 
power/energy  limiters.  Defocuslng  media  such  as  Si,  InSb  and  CdxHg1_xTe  are 
generally  preferable  to  focusing  media,  since  in  the  latter  cases  the  focal 
spot  will  rapidly  move  towards  the  focusing  system,  leading  to  a  possibility 
of  damage  to  the  device  rather  than  the  reverse.  some  of  the  advantages  of 
this  type  of  limiting  device,  when  compared  with  other  methods  of  limiting, 
are  (l)  complete  passivity,  (2)  fast  response  and  recovery  times,  (3) 
broadband  operation,  (4)  absence  of  low-power  optical  filtering,  (5)  large 
dynamic  range.  Broadband  operation  i3  useful  where  the  incoming  radiation  is 
of  uncertain  wavelength  or  where  the  laser  is  frequency-agile.  Remaining 
problems  which  should  be  addressed  in  the  future  include 

(1)  Materials  research  directed  to  obtaining  large  values  of  n2  and  o^, 

(2)  Production  of  high-quality  media  for  imaging  applications, 

! 3 )  Study  of  the  effects  of  coherence  quality  of  the  beam  (or  pulses), 

( 4 )  Methods  for  minimising  possible  turbulence  at  low  powers. 


Lastly,  the  recent  experimental  and  theoretical  research  by  Khoo  et 
al  (18)  and  Khoo  (19)  on  nonlinear  refraction  of  a  gaussian  laser  beam  by  a 
thin  nematic  liquid-crystal  film  has  demonstrated  both  ringed  structures  in 
the  far-field  intensity  distribution  and  strongly  power-dependent  oscillations 
along  the  optic  axis.  These  features  are  in  agreement  with  the  predictions 
of  our  model,  and  lend  support  to  the  choice  of  refractive  nonlinearity 
(quadratic  in  terms  of  the  susceptibility,  or  cubic  in  terms  of  the 
macroscopic  polarisation)  where  liquid  crystal  films  are  employed  as  the 


a 


nonlinear  medium.  Essential  procedures  for  incorporating  additional 
I  focusing,  lenses,  and  mirrors  into  the  model  have  been  discussed  in  detail  by 

these  authors.  The  conversion  from  the  model  represented  in  fig.  la  to  that 
shown  in  fig.  lb  is  relatively  straightforward.  Related  work  with  self-phase 
modulation  by  nematic  films  also  has  been  carried  out  very  recently  by 
r  Santamato  and  Shen  1201.  However,  it  should  be  noted  that  despite  their  very 

>  large  refractive  nonlinearities,  the  present  liquid-crystal  systems  appear  to 

be  far  too  slow  to  be  considered  for  the  purpose  of  limiting  pulses  of 
r  nanosecond  or  shorter  duration. 

1 
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-  - *  Defocusing  IS 

- Focusing 

— - -  Low  Power 


(a) 


(b) 


Arrangement  for  a  self-focusing/defocusing  power  limiter: 

(a)  simple  model,  (b)  practical  model.  The  nonlinear  medium  is 
designated  NL,  and  the  screen  s. 


FIGTOE  2  Normalised  intensity  distribution  for  the  gaussian  input  F(x)  - 
exp(-x2)  in  the  region  of  the  linear  focal  plane.  Radiation  is 
incident  from  the  left  for  self-focusing  media  and  from  the  right 
for  self-defocusing  media.  Fig.  (d)  is  a  close-up  of  Fig.  (c) 
with  an  expanded  scale. 


Focal  Plane 


FIGURE  2(d) 


Iaophotes  near  the  linear  focal  plane,  for  gaussian  input 
amplitude;  (a)  P  -  o,  (b)  p  -  5.  Intensity  maxima  are  designated 
H  and  minima  L.  contour  Intervals  are  half  orders  of  magnitude; 
the  highest  contour  In  (a)  is  10-1  2  and  In  (b)  is  10*1. 
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FIGURE  5(b) 


6  isophotes  near  the  linear  focal  plane,  for  parabolic  Input 
amplitude;  (a)  p  -  o,  (b)  p»  5.0,  (c)  p  -  10.  Maxima  are 
designated  H  and  minima  L.  contour  intervals  are  spaced  on  a 
semi-log  scale  with  ten  equal  divisions  per  order  of  magnitude. 
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FIGTOE  8  Transverse  profile  for  gaussian  input  at  p  -  8.02:  (a)  Full  line, 

experimental  (ref.  l);  (b)  Broken  line,  theoretical  matched  at 

t  -  0.6. 
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